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Abstract 

We discuss photon Bremsstrahlung induced by virtual graviton exchange in proton-proton 
interactions at hadronic colliders, resulting from the exchange of Kaluza-Klein excitations 
of the graviton. The relevant subprocesses, gg — > G — > e + e"7 and qq — > e + e~7 are 
discussed in both the ADD and the RS scenarios. Although two-body final states (or real 
graviton emission) would presumably be the main discovery channels, a search for three- 
body final states could be worthwhile since such events have characteristic features that 
could provide additional confirmation. In particular, the /c^-distribution of the photon is 
in both scenarios harder than that of the Standard-model background. 
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1 Introduction 



The idea of additional compact dimensions and strings at the TeV scale, proposed by 
Antoniadis pQ for solving the hierarchy problem, together with the idea that Standard- 
model (SM) fields live on branes in a higher-dimensional space j2] have led to the even 
more radical speculations that extra dimensions might be macroscopic, with SM fields 
confined to the familiar four-dimensional world (brane) [3J H] . The models which allow 
for gravity effects at the TeV scale can be grouped into two kinds: those of factorizable 
geometry, where the extra dimensions are macroscopic [Sj ("ADD scenario"), and those 
of non-factorizable (warped) geometry, the simplest example of which has only one extra 
dimension separating "our" brane from a hidden brane [3] ("RS scenario"). 

In both these scenarios, the propagation of gravitons in the extra dimensions leads to 
gravitons which from the four-dimensional point of view are massive. In the ADD scenario, 
these Kaluza-Klein (KK) gravitons have masses starting at values of the order of milli-eV, 
and there is practically a continuum of them, up to some cut-off Ms (close to the effective 
Planck scale) of the order of TeV, whereas in the RS scenario they are widely separated 
resonances with mass splittings of the order of TeV. In both cases, they have a universal 
coupling to matter and photons via the energy-momentum tensor. 

These recent speculations have led to several studies |HHT3] of various experimental 
signals induced by graviton production and exchange. The new scenarios allow for the 
emission of massive gravitons [SI E] ; which would lead to events with missing energy 
(or transverse momentum), as well as effects due to the exchange of virtual gravitons (in 
addition to photons and Z's) El IH1 AT] . These processes include the production of 
dileptons and diphotons in electron-positron collisions, as well as gluon-gluon and quark- 
antiquark-induced processes at the Tevatron and LHC. 

In fact, several searches at LEP and the Tevatron have given direct bounds on the 
effective Planck scale, of the order of a TeV [T3Tll5| . while astrophysical arguments result 
in very strong limits when applied to the simplest ADD scenarios, for n = 2 and 3 extra 
dimensions [IS]. Of course, the direct experimental searches are most worthwhile. The 
above studies all focus on two-body final states, which are expected to be dominant, and 
therefore lead to the most stringent bounds on the existence of extra dimensions. 

Here, we shall investigate photon Bremsstrahlung induced by graviton exchange [T7] . 
While this cross section is further reduced by 0(a/7c), so is the background. It has some 
characteristic features resulting from the exchange of a spin-2 particle and from the direct 
graviton-photon coupling, that we would like to point out. These features may be useful 
in discriminating any signal against the background. 

Specifically, we shall consider the process 

pp — > e + e~7 + X, (1.1) 

which may get a contribution due to graviton exchange, and which for energetic electrons 
(or muons) and photons should experimentally be a very clean signal. (There is also a 
related process, where a graviton is emitted in the final state [T8].) 
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Since this final state is very distinct, and since the Standard-model (Drell-Yan) back- 
ground is well understood, the process (jl.lj) may offer some hope for observing a signal or 
improving on the exclusion bounds. 

This paper is organized as follows: First (Sec. EJ) we consider the gluon-gluon fusion 
contribution to both two-body and three-body final states. Then (Sec. EJ) we consider 
quark-ant iquark annihilation, which also gives rise to the Standard-model background. 
We calculate the cross section by summing over the KK tower within the ADD (Sec. EJ 
and RS (Sec. EJ) scenarios for a selected choice of parameters, and finally we give some 
concluding remarks (Sec. El). 

2 Gluon— gluon fusion 

We shall first discuss gluon-gluon fusion. Due to increasing gluon luminosity at high 
energies (LHC), this contribution will be dominant for a certain range in invariant mass 
of the (three-body) final state. At even higher invariant masses, the contribution to the 
signal, from quark-antiquark annihilation, becomes dominant. 

2.1 Two-body final states 

The process of interest, Eq. (jl.lj) . is related to the two-body final state 

pp -> e + e- + X, (2.1) 

which may proceed via gluon-gluon fusion and an intermediate graviton, 

gg - G - e+e". (2.2) 

For massless electrons, the cross section for single graviton exchange resulting from 
gluon-gluon fusion is (in agreement with the results of 



U 99^ee ^ _ + ( m .p.)2' ^ '°) 

with s = (fci + /c2) 2 the two-gluon invariant mass squared. Furthermore, and are the 
mass and width of the graviton 1 , and k is the graviton coupling, to be defined below. The 
angular distribution, which is forward-backward symmetric, is given by 1 — cos 4 #, where 
6 is the cm. scattering angle. 

With £i and £2 the fractional momenta of the two gluons, k\ = £i-Pi, ^2 = £2-^2, and P\ 
and P2 the proton momenta, (Pi + P2) 2 = s, we have s ~ For the over-all process 

(12. lj) we thus find the differential cross section 

da (G) f 1 f 1 da {G) 



Jo 



ds J q Jo ds 

1 The Kaluza-Klein index, n on and should not be confused with n, the number of extra 
dimensions. 
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= / dfl / dfc /,(&)/,(&) ^ (66* - 5) <^-lee(*) = ~^(*) ^ee(S), (2.4) 

Jo Jo s 

with the relevant convolution integral, I gg (s), over the gluon distribution functions given 
by Eq. (|A.1J) in Appendix A. 

2.2 Three-body final states 

Let us now consider the contribution from gluon-gluon fusion to the Bremsstrahlung pro- 
cess in Eq. (jl.lj) . The underlying subprocess, 

gg _> G _ e + e - 7) (2.5) 

can proceed via the four Feynman diagrams of Fig. ^ the basic couplings for which are 
given by Han et al. [2j (see also Giudice et al. jHj). 




Figure 1: Feynman diagrams for gg — > G —>■ e + e~ r y. (There is a corresponding set of 
diagrams with qq initial states, which we shall refer to as "set B" . The SM diagrams 
corresponding to (1) and (2), with qq initial states and a photon or Z instead of the 
graviton, will be referred to as "set _D 7 " and "set Dz" •) 

The evaluation of the cross section is straightforward, and the differential cross section 
(w.r.t. the azimuthal angle, \i an d cos#, where 9 is the angle between the photon and the 
beam) is of fourth order in the invariants. This is due to the underlying mechanism being 
the exchange of a spin-2 object. It is straight-forward to verify that it is gauge invariant 
with respect to the gluons, as well as to the photon (actually, diagram 4 is by itself gauge 
invariant). But the expression is quite lengthy, so we shall not write it out here. 
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The angular distribution of the (two-body) non-radiative cross section ()2.3|) is given by 
fourth-order polynomials in cos#. Here, just like in gluon Bremsstrahlung (see, e.g., [TT?]). 
there is an accompanying dependence on the azimuthal angle Xi but now up to fourth 
order in cosx, or, equivalently, up to cos 4%. 

After integrating over the azimuthal orientation of the events, the cross section is of 
the form 



gg^ee-y 



, , , ,, b (x 1 ,x 2 ) + b 2 (x 1 ,x 2 ) cos 2 9 + bJx 1 ,x 2 ) cos 4 6, (2.6) 
dx\dx 2 d cos6 

similar to the two-body final states, i.e., the gluon-gluon fusion does not contribute to any 
forward-backward asymmetry. 

In our calculations we have chosen the unitary gauge (C,^ 1 = in the notation of 
[7j), whereby the scalar exchange decouples. After averaging and summing over gluon, 
electron and photon polarizations, and integrating over event orientations w.r.t. the gluon 
momentum, we find (for exchange of a single graviton) 

d 2 rr {G) a^nO 2 s 2 

dx x dx 2 2560tt 2 (s - m 2 ) 2 + (m H V n ) 2 B{ 1: 2h K ] 

with Xb{x\,x 2 ) given by Eq. (jA.8|) in Appendix A. Furthermore, a is the fine-structure 
constant, Q e = —1 is the electron charge, and x\, x 2 and X3 denote the fractional energies 
of the electrons and the photon in the cm. frame, 

xi = E\/V§, x 2 = E 2 /V§, X3—U/VI, < Xi < |, (2.8) 

with x\ +x 2 +X3 — 1. The expression Xb{xi, x 2 ) exhibits the familiar singularities in the 
infrared and collinear limits, s\ = (pi + k) 2 = s(l — 2x 2 ) —>■ 0, s 2 = (p 2 + k) 2 = s(l — 2xi) — > 
0, as well as a collinear singularity at s 3 = {p\ + p 2 ) 2 = s(l — 2x 3 ) — > due to the fourth 
Feynman diagram. Here s = (ki + k 2 ) 2 = (pi +p 2 + k) 2 . The additional singularity means 
that there is a tendency to have events with hard photons J7j. This is one way in which 
these events differ from ordinary QED-based Bremsstrahlung. 

The differential cross section in Eq. ()2.7j) can be written more compactly as 

1 d 2 a { g % ei _4aQ 2 e (2 g) 



n-( G ) dX\dx 2 7T 

with agg-^ee given by Eq. (j2.3J) . As we see, the cross section is reduced by a factor (D(a/7r) 
compared to the two-body cross section. 

Analogous to Eq. ()2.4|) . we find for the gluon contribution to the over-all process (jl.lj) 

U Ogg^ee^y _ ±j » gg^ee-y ,^ ^qn 

dsdx\dx 2 s 99 dx\dx 2 ' 
with the convolution integral given by Eq. (|A.1J) in Appendix A. 
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3 Quark— ant iquark annihilation 



Another process which contributes to (jl.lj) . and gives rise to the SM background, is quark- 
antiquark annihilation. Furthermore, at large invariant masses of the final state, it also 
gives the most important contribution to the signal. 



3.1 Two-body final states 

The process in Eq. (12. 1J) may also proceed via quark-antiquark annihilation and an inter- 
mediate graviton, with the following cross section for single graviton exchange (initial state 
quarks are considered massless) 

a (o) = g2 f3 d 

°m->ee 1536Q7r (s_ m 2)2 + m 2 r 2> ^ 

in agreement with Ej- It differs from the cross section for gluon-gluon fusion by a 
factor 2/3. 

There is also a SM background to this process, where the same final state is produced 
through photon or Z exchange. This well-known cross section is given by 

a^tl(s) = [Q 2 q Ql + 2Q q Q e v q v e Re X (s) + {v\ + a 2 q )(v 2 e + a 2 e )\ X (s)\ 2 ] , (3.2) 

with 

1 § 

X{S) = sin 2 (2M (s - m\) + im z T z ' (3 " 3) 

where we have normalized vector and axial- vector couplings to Vf = Tf — 2Q f sin 2 9w and 
cif = Tf respectively, with Tf the isospin. Furthermore, mz and V z are the mass and width 
of the Z boson, Q q the quark charge, and 9w the weak mixing angle. 

In the case of qq — > G — > e + e _ , with the cross section given by Eq. (|3.1|) . the angular 
distribution is 1 — 3 cos 2 6+ A cos 4 9, whereas for the photon exchange, Eq. (|3.2jl . the angular 
distribution is given by the familiar 1 + cos 2 9. The interference between graviton and 
photon exchange has an angular distribution given by cos 3 9 (as pointed out by Ref. i.e., 
it exhibits a forward-backward asymmetry and vanishes upon integration. The interference 
between graviton and Z exchange exhibits a slightly different angular distribution (which 
also vanishes upon integration). 

For pp collisions, we find the graviton contribution to the differential cross section (in 
accordance with Eq. (|2.4jl ) 

dn {G) 1 

-^^ = - s I B B(s)a^l ee (s), (3.4) 

with Ibb{s) given in Appendix A. The SM contribution can be found in a similar manner, 
but here the convolution integrals must be weighted by the factors Q q , Q q v q and (v 2 + a 2 ) 
for photon exchange, interference between the photon and the Z, and for Z exchange, 
respectively (see Appendix A). The reason for this is that the convolution integral implicitly 
contains flavor summation. 
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3.2 Three-body final states 



Now we will examine the subprocess 

qq — > e + e~7, (3.5) 

which is determined by four sets of Feynman diagrams. First, there are the diagrams 
of Fig. El referred to as "set A" , describing "initial-state radiation" . Then, there are four 
diagrams analogous to those of Fig. [TJ where the initial-state gluons are replaced by quarks 
and antiquarks. We shall refer to these as "set B" , they describe "final-state radiation". 
Finally, we have the SM background, which arises from diagrams similar to 1 and 2 of 
sets A and B, with a 7 or Z exchanged instead of the graviton. We shall refer to the 
SM diagrams as "sets C 7 ", "Cz" (initial-state radiation), "-D 7 ", and u Dz" (final-state 
radiation). It is convenient to separate initial- from final-state radiation, since, in the 
former case, the propagator does not carry all the momentum of the initial quarks. 




Figure 2: Feynman diagrams for "initial-state radiation" in qq — > e + e~ r ). We refer to these 
diagrams as "set A" . The SM diagrams corresponding to (1) and (2), but with a photon 
or Z instead of the graviton, are referred to as "set C 7 " and "set Cz" . 



The gg-initiated cross section can be decomposed as 

-(G) .(SM) ~(G, 7 ) , ~(G,Z) f o fi N 

u qq^ee'y > u qq— *ee-y > u qq^ee^f > emu. <J qq^ ee ^ , ^<j.uj 

where the first term is the graviton contribution (sets A and B), the second term is the 
Standard-model background (sets C and D) and the last two are graviton-photon and 
graviton-Z interference terms, respectively. 
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frifLeej = &A + &B + &AB- (3.7) 



First we shall consider the graviton exchange diagrams, where we introduce the follow- 
ing notation, 

^ qq— >ee7 

For initial-state radiation (set A), we find 
d 2 a A an 4 s 3 Q 2 



dxidx 2 12288tt 2 (s 3 - ml) 2 + (m H T f 



X A (xi,x 2 ), (3A 



where X A (xi,x 2 ) is given by Eq. (|A.4|) in Appendix A. Since the denominator in Eq. ()3.8j) 
depends on s 3 = (1 — 2x 3 )s instead of s, this contribution will be smeared out when 
integrated over x 3 . 

The corresponding result for final-state radiation (set B), is 

d 2 a B a^sQ 2 e s 2 

X B {xi,x 2 ), (3.9) 



dx x dx 2 3840vr 2 (s - m|) 2 + {m H T n ) 

with X B (xi, x 2 ) given by Eq. (jA.8|) . This contribution is, like in the two-body case, identical 
to the contribution from gluon-gluon fusion, except for a factor of 2/3. 

There is also an interference term, &ab, between sets A and B, which has to be con- 
sidered. It contributes to the forward-backward asymmetry, and vanishes when integrated 
over all event orientations. 

In the three-body case, the SM background becomes 

crfq^ eei = &c + &D + &CD- (3.10) 

where the contribution of initial- and final-state radiation is given by 

d 2 &c 3aQ 2 ^ (SM) 
dx^ 2 =— a ^ ee{S3)XciXuX2) > 



d 2 &D ^aQl A 



dX\dx 2 7T 



Ke(«) X D {x Xl X 2 ), (3.11) 



with Xc{xi,x 2 ) and Xd(xi,x 2 ) given in Appendix A, and Ogf^ee given by Eq. (|3.2)1 . This 
<td is the familiar Bremsstrahlung cross section expressed in terms of the related two-body 
process. While the contribution of the photon-exchange part of the interference between 
initial- and final-state radiation to the integrated cross section, cfq 1 d 1 i vanishes [20J, this 
is not the case for terms involving Z-exchange |21| . They are included in the full SM 
background, Eq. (J4.8|) . 

For the interference terms between graviton exchange and the SM diagrams, we intro- 
duce the following notation: 

~(G, 7 ) ~ ~ ~ ~ 

<7 q q^ en — a AC-y + & AD 1 + &BC 1 + ®BD 1 - I 

(J q qL, eei = <J ACz + cr ADz + a BCz + Obd Z i (3-12) 
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where the subscripts indicate the diagram sets involved. 

We find that both ctac 1 an d <5"rd 7 , together with the terms of &ac z an d &bd z which are 
proportional to the vector coupling, vanish after integration over event orientations, but 
they contribute to the forward-backward asymmetry. The oag z term proportional to the 
axial coupling also vanishes upon integration over all event configurations. Finally, after 
integration over angles the obd z term proportional to the axial coupling, being propor- 
tional to the Legendre polynomial P2(cos8), also vanishes. (When introducing cuts, the 
cancellation is not complete.) 

For the non-vanishing graviton-SM interference terms, we get 



d 2 a A D-, 
dx\dx 2 

d 2 OAD z 

dx\dx 2 

d 2 VBC\ 

dx\dx2 
d 2 o BCz 
dx\dx 2 



o?K 2 Q 2 q Q\ ^ 

a 2 K 2 Q q Q e v q v t 

144vr 
a 2 K 2 Q 2 Q 2 



s 3 



s 3 -mt 



Re 



s 3 



S3 



mt + vm 



X AD (x 1 ,x 2 ), 

X AD {x 1 ,x 2 ) 



Re 



1447T 

a 2 K 2 Q q Q e v q v, 



1447T 



s -ml + irrifir H 



Xbc(xi,x 2 ) 



X*(s 3 ) 



s — ml + im H T r 



X BC (x 1 ,x 2 ) 



(3.13) 



where x(s) is given by Eq. ()3.3|) . whereas Xad(xi,x 2 ) and X BC (xi,x 2 ) are given in Ap- 
pendix A. 

To find the contribution from quark-ant iquark annihilation to the over-all process (jl.lj) . 
a relation similar to the one given by Eq. ()3.4j) should be used. However, when there are 
quark charges or vector /axial- vector couplings (v q or a q ) involved, the convolution integrals 
must be weighted by these factors, as shown in Appendix A. 



4 Bremsstrahlung in the ADD scenario 

In this section we shall consider Bremsstrahlung within the ADD scenario [3]. First, we 
shall present the cross section as a function of invariant mass, and then study the photon 
distribution (or k±_ spectrum) of such events. 

4.1 Total cross section 

In the ADD scenario, the coupling of each KK mode to matter is Planck-scale suppressed. 
However, since the states are very closely spaced, with [2j 

< - ^f, (4.1) 

and -R/2-7T the compactification radii, the coherent summation over the many modes leads 
to effective couplings with strength I /Ms, where Ms is the UV cut-off (close to the effective 
Planck scale). 
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Explicitly, in this scenario, the graviton coupling is in the (4 + n)-dimensional theory 
given by [7j 

9mn = Vmn + kh MN , k 2 = 167rG^ 4+n) , (4.2) 

where G^ +n ^ is Newton's constant in 4 + n dimensions. In 4 dimensions the coupling can 
be written as 

k 2 = V^k 2 = WnV- 1 G ( * +n) = 16ttG n , (4.3) 

with V n the volume of the n-dimensional compactified space (V n = R n for a torus T n ) and 
Gn the 4-dimensional Newton constant. 

Summing coherently over all KK modes in a tower, the graviton propagator gets re- 
placed 



1 



s -ml + imnTf. 



-%D{s) 



R n 



r(n/2) (4tt)"/ 2 



[2/(M £ 



VK\ 



(4.4) 



with 



I(M< 



n/2-l 



V — ( ^ 

^ 2k U/a 



fc=i 

(n-l)/2 

E 

fc=i 



1 



2A;- 1 



2fc 



Ms 



. Mo 
-log -^-1 , 



2fc-l 



M 5 + 



Mo- 



ri = even, 



n = odd. 



(4.5) 



for n extra dimensions. 

Higher order loop effects may be important [22], so these expressions should not be 
taken too literally. In particular, this applies to the dependence on the number of extra 
dimensions. In the approach of [Sj and [Hj this uncertainty, including the n-dependence, is 
absorbed in the cut-off in such a way that D(s) and D(ss) (see Eq. ()4.4j) ) are indistinguish- 
able. Here, in order to preserve the qualitative difference between these two propagators, 
related to final- and initial-state radiation, we shall use the expressions of Eq. (j4.5J) . 

Invoking the relation [7j 



2 nn 



K Z R 



^(Aix) n / 2 T(n/2)M~ {n+2) 



(4.6) 



between the volume of the extra dimensions, the gravitational coupling and the cut-off 
scale, the differential cross section can be expressed as 



d 3 a 



as 3 



dsdx\dx2 192sM| 
olQ\s 



I DjDj (s) [4I 2 (M s /^) + ir 2 ] X A (x h x 2 ) 



+ 



120sMi \M\ 



-2 [3I gg (s)+2I BB (s)\ AI 2 (M S /Vs) + 



IT 



X B (xi,x 2 ) 



a 



18sM 2 L 



n/2 

I ( 2( 3eW> 



{§) I(M S /yfe 
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+ QeV e I Dl D z (s) Re{ X *{s) [2/(M 5 /Vii) - wr] })X AD {x 1 ,x 2 ) 

(a \ n/2 
(2Q2j D7D7 ( S )j(M 5 /v / I) 



+ Q^Jd,d z {s) Re{x*(s 3 ) [2I(M S /V§) - m] }^X BC (x u x 2 ) 

+ dV SM > 
ds dx\dx 2 ' 

where x(§) * s given by Eq. ()3.3j) and the convolution integrals are given in Appendix A. 
The SM background is given by 

dV SM) a 3 

du^r = W^xM 

+ (^e + a e) / C z C z (s)|x(s 3 )| 2 } A C (xi,X 2 ) 

+ -^{QIId^(§) + 2Q e vJ DjDz (s) Re X (§" 



+ {v 2 e + a 2 e )I DzDz {s)\x{s)\ 2 } X D (x u x 2 ) 
— ^—-^{Qe Ic^d z {s) Rex(s) + QelczD^s) Rex(ss) 

oSS 

+ Av e I Cz D z (s) Re[x*(s)x(s 3 )}} X C d(xi, x 2 ). 



(4.1 



In Eq. (|4.7jl . the different contributions are given in the following order: First quark- 
antiquark annihilation with initial-state radiation (set A), then gluon-gluon fusion and 
quark-ant iquark annihilation with final-state radiation (set B), various interference terms 
between graviton exchange and SM amplitudes, and finally the SM background. The 
origins of the different terms are reflected in the subscripts of the Vs. 

In Fig. El we show the cross section, differential w.r.t. VI (labeled M ee7 in the figures), 
for Ms = 4 TeV and n = 2,3,4,5 and 6 (left panel) in the ADD scenario, where we 
have integrated out the Xi,x 2 dependence (see Appendix B). The right panel shows only 
the n = 2 curve, with the contributions from gluon-gluon fusion and quark-antiquark 
annihilation (induced by graviton exchange) also displayed. Note that gluon-gluon fusion 
is dominant from vl ~ 1.3 TeV up to Vs ~ 3 TeV for this choice of parameters, whereas 
the quark-antiquark annihilation process takes over at larger invariant masses. Below 
~ 1 TeV, the background is larger than the signal. 

We have integrated over x™ in < x 3 < 0.5, subject to the y-cuts: Si,s 2 > ys, s 3 > 
y 3 s, where both y = y 3 = 0.01. The minimum invariant mass of the two electrons is 
controlled by y 3 . At a scale \/j = 1 TeV, the cut of y = 0.01 corresponds to electron (or 
photon) energies exceeding 10 GeV. We consider a minimum x 3 of 0.1 in these plots. The 
corresponding resolution is well within that foreseen at the LHC j2B] ■ For the angular cuts, 
we take the pseudorapidity \r]\ < 2.5. 

Here, and in the remaining figures, we use y/s = 14 TeV, which corresponds to the 
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Figure 3: Cross sections for pp — > e + e~7 + X at y/s = 14 TeV. Both plots are for the 
ADD model with Ms = 4 TeV. We have set the number of extra dimensions (from above) 
to n — 2, 3, 4, 5, and 6 for the left panel. The full differential cross section, da/d^/I (solid), 
and the SM background (dashed) are shown in both plots, whereas the gg (dash-dotted) 
and the qq (dotted) contributions to the signal are shown in the right panel for n = 2. 



LHC energy. With an integrated luminosity of 100 fb~ and a bin-width of 100 GeV, we 
might expect a few events per bin at invariant masses above 1 TeV. 

Near the Z mass, the cross section has an additional peak (barely visible in the figures) 
since both the amplitudes with s (set Cz) and those with S3 (set Dz) resonate. Also, at 
s ~ 10 x rriz ~ 0.9 TeV, the cross section has a fall-off. In fact, both the second peak 
and the fall-off are related to "radiative return", where the emitted photon has the right 
energy to make the Z propagator resonate. The peak near mz is at the starting point for 
radiative return, determined by S3 (and hence the lower value for £3), whereas the fall-off 
near 0.9 TeV is related to the end point of the radiative return, given by the upper value 
of X3, which is ~(1 — y| ut ). 

Close to the cut-off, Ms, the cross section blows up due to the logarithm in I (Ms /VI). 
This is of course an artifact, due to the way the cut-off is treated [22J. Note that the 
explicit n-dependence in Fig. El is in the approach of 0|H] absorbed in the cut-off. 

Let us now discuss the different terms related to qq annihilation. The part of the qq 
annihilation cross section which is related to Xa (initial-state radiation with graviton ex- 
change) is everywhere significantly smaller than the one related to Xb (final-state radiation 
with graviton exchange), by more than an order of magnitude. This is partly due to the 
difference in the convolution integrals. 

Among the different interference terms, related to Xad and Xbc, the most impor- 
tant one is that between set B (final-state radiation with graviton exchange) and set C 7 
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(initial-state radiation with photon exchange). For invariant masses below ~ 1 TeV, this 
interference dominates the qq part of the signal (but here, it is overwhelmed by the SM 
background). We note that the SM amplitude does not interfere with the gluon-gluon 
fusion part of the graviton-mediated amplitude. 

4.2 Photon distribution 

Whereas photons emitted by QED Bremsstrahlung tend to be soft and/or collinear with 
the fermions, the present graviton-induced Bremssstrahlung will be more energetic [T7j . 
and also emitted at larger angles. A first, qualitative manifestation of this feature is that 
the transverse momentum spectrum (with respect to the incident beam) will be harder. 
We shall here make this statement quantitative. 

In Sec. 14.11 we have integrated the differential cross section over Xi,x 2 and cos 9 to 
obtain cross sections which are differential w.r.t. y/§ (see also Eq. (|2.4jl ). Now we will 
instead change variables from (xi, x^, cos#) via (xi% = %i — X2,xs,cos9) to (xu, k±, k\\), 
where k± = usmO and k\\ = ucos9. Here 

+ fcjj = k = u = x 3 V§, (4.9) 

with uj the energy of the photon. After integration over £12, k\\ and y/§, we get the cross 
section, differential w.r.t. k±, the photon momentum perpendicular to the beam. 

In Figure left panel, we show da/dk± for different numbers of extra dimensions, and 
with Ms = 4 TeV, where we have integrated over Vs up to the cut-off, Ms- In accordance 
with Fig. El we impose the cuts: x| ut vl < k < |(1 — 1/3) Vs. Since the photon should not 
be too close to the beam axis, we also require k± > 100 GeV. We see that the fc^-spectrum 
becomes much harder when extra dimensions are involved. 

The right panel of Fig. H] shows bin-integrated cross sections. Here we have chosen 
invariant masses of \f§ = 1,2,3 TeV, and integrated d 2 a/(dk±d\/I) over 100 GeV bins 
around these values. The peaks at 0.1, 0.2 and 0.3 TeV are also related to cuts on k (or 
23). Note that the SM background is very small at higher invariant masses. 

5 Bremsstrahlung in the RS scenario 

As in the previous section we will also for the Randall-Sundrum scenario [1] first present 
the total cross section, and thereafter the photon distribution. 

5.1 Total cross section 

In the RS scenario, the graviton masses are given by |TU] 

m n = kx n e~ kr ^ = — mi, (5.1) 

Xi 
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1.5 2 0.5 1 1.. 

kJTeV] kJTeV] 



Figure 4: Photon distributions for the ADD model. Left panel: da/dk± (solid) with 
(from above) n = 2,3,4,5,6 extra dimensions and Ms = 4 TeV, together with the SM 
background. Right panel: da/dk± (solid) for n = 2 and Ms = 4 TeV, integrated over 
100 GeV bins around M ee7 = M = 1, 2, 3 TeV, with the corresponding SM backgrounds 
(dashed). 



where x n are roots of the Bessel function 2 of order 1, Ji(x n ) = 0, k is of the order of the 
(four- dimensional) Planck scale and r c the compactification radius of the extra dimension 3 . 
Since there is only one extra dimension in this scenario, we shall use m n instead of m^ for 
the mass of the nth graviton. 

The gravitational coupling is in this model given by fUl ^] 

k = y/lfa — = V2 ^- J-, M P1 = ^ = 2.4x 10 18 GeV. (5.2) 



m n Mpi mi Mpi 



57T 



Here we shall use mi and k/Mp\ as the two parameters which specify the model. Note 
that 0.01 < fc/Mpi < 1 BUj. 

The differential cross section can in the RS scenario be expressed as (with the different 
contributions given in the same order as in Eq. ()4.7jl ): 

d 3 cr as 3 ( x x \ ( k \ 2 , 

Id~d~{s) \G{s 3 )\ X a (x 1 ,x 2 ) 



dsdx\dx2 30727r 2 s \ m i/ \M 



■ pi 



+ A fe) 4 (wS Q ' l3I " iS) + 2/BB(i)l |5(i)|2 Xb{xuX2) 



2 The first four roots are 3.83, 7.02, 10.17 and 13.32. 
3 To solve the hierarchy problem, kr c ~ 12 is required. 
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a / x\ \ ( k 



72ns \mi J \^pi 

X 



QllD^is) Re£(s 3 ) + Q e v e I DiDz (s) Re [x*(s) Q(s 3 )} JX AD (x 1 ,x 2 ) 
+ (Q 2 e I D ^(s)Reg(s) +Q e vJ DiDz {s) Re [x*{s z )Q{s)\ )x BO (x 1 , x 2 ) 

d 3 a (SM) 

ds dx\dx2 ' 



(5.3) 



where the SM contribution is given by Eq. (|4.8jl . For massive graviton exchange we have 
introduced the abbreviation 



^s-mH im n T, 



with 121 



Tn = m < * = m^{i&) ' (5 - 5) 



and 



7g = 1 + X 7 + Xz + Xw + Xe + Xq + Xh + Xr (5.6) 
the total graviton width in units of the two-gluon width. Neglecting mass effects, we have 

1 13 13 



Xl 8' XZ 96' XW 48 



N e N c N q 1 

Xe = —, X q = —^r-, XH = Xr = ^- (5.7) 

Here, = 6 is the number of leptons, and N c N q = 18 is the number of quarks weighted 
with color factors. Note that since we have neglected quark and electron masses, there 
is no contribution to the cross section from radion exchange, since the radion couples to 
the trace of the energy-momentum tensor. However, it contributes slightly to T n through 

Eq. iprsji . 

We display the RS scenario cross section, differential w.r.t. \/~§ (see Eq. ()5.3|) ). in Fig.0 
In the left panel, we have summed over KK states, and chosen the first graviton resonance 
at rrii = 1 TeV, with fc/Mpi set to 0.1, 0.05 and 0.01 (from above). In the right panel, 
we show the different contributions (gg and qq) to the cross section (for k/Mp\ = 0.05) 
induced by graviton exchange. The cuts are the same as in the ADD case. 

As we mentioned above, the integration over S3 (or £3) will smear out the contribution 
from initial-state radiation. For invariant masses greater than 

mi 

1.1 mi, (5.8) 



'1 - 2xf m 



s™ ax will always be greater than mi, and radiative return to the first graviton resonance 
gives a small peak at this value, ()5.8j) . and similarly for the higher resonances. In fact, these 
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Figure 5: Differential cross sections for the RS model, with m\ = 1 TeV. The left panel 
shows the differential cross section, dajd\fk (solid) with (from above) k/Mp\ = 0.1, 0.05 
and 0.01, and the SM background (dashed). Right panel: the gg (dash dotted) and qq 
(dotted) contributions to the cross section are shown for k/Mpi = 0.05. 



secondary peaks are mainly due to the term involving Xa (initial-state bremsstrahlung with 
graviton exchange), which gives no visible effect in the ADD case. 

Similar to the ADD case, for invariant masses below ~ 1 TeV, it is the interference 
between set B and set C 7 which dominates the qq signal (dotted). Between the peak at 
1 TeV and the secondary peak next to it (due to radiative return) this interference term 
becomes negative, resulting in a dip in the total cross section. In fact, the total qq signal 
(graviton exchange plus its interference with the SM amplitude) is negative in some small 
region of invariant mass. However, these structures depend on details of the cuts imposed. 

It should be noted that the third and fourth peaks should be somewhat reduced since 
we have not taken into account that these gravitons can decay into the first KK resonance. 
Self- interactions of the gravitons were considered in where a BR of about 15% was 
found for the G 3 — > G\G\ decay. 

By comparing Figs. El and we see that the difference between the two scenarios is 
striking, with extremely narrow, and widely separated resonances in the RS case, compared 
to the continuum in the ADD case. Note that according to our expressions, the k/Mp\ 
dependence cancels at the resonance. A single peak should therefore have a height which is 
independent of k/M P y, but at very high invariant masses we see that the fc/Mpj = 0.1 peak 
is higher than the other peaks, and also slightly shifted. This is mainly due to interference 
with, and overlap of the neighboring peaks which are very broad. 

Since the cross sections for graviton-induced Bremsstrahlung are much lower than for 
the corresponding two-body final state, an important question is, however, if there is any 
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Figure 6: Integrated cross sections for the RS scenario. In the upper left panel we have 
integrated the cross sections in Fig. |5]over for > M™" 1 . The choice of parameters 
is k/Mpi = 0.1, 0.05 and 0.01 (from above) for m x = 0.5 (solid), 1 (dotted) and 1.5 TeV 
(dash-dotted). The SM background (dashed) is also shown. The other three panels show 
the corresponding result of integrating over 100 GeV bins, with ni\ = 0.5 (upper right), 
rrii = 1 (lower left) and mi = 1.5 TeV (lower right). 



chance of seeing these resonances in the experiments. To give an order of magnitude 
estimate of the number of events to expect from these narrow peaks, we have integrated 
the differential cross sections given in Fig. El over bins in M en . In the upper left panel of 
Fig. IHlwe integrate over \f§, starting from = M™^. The different curves correspond 
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to fc/Mpi = 0.1, 0.05 and 0.01 (from above) for m x = 0.5 (solid), 1 (dotted) and 1.5 TeV 
(dash-dotted). This should be compared to the SM background (dashed) which is also 
shown. In the remaining panels of Fig. El we have integrated over 100 GeV bins for the 
same choice of k/Mp\, for mi = 0.5, 1.0 and 1.5 TeV. 

We see that for these parameters it should be possible to resolve at least the first peak, 
and in most cases several peaks are visible. We emphasize that this is not a simulation, 
but these plots should provide an indication of the number of events that correspond to 
these narrow peaks. 



5.2 Photon distribution 

The photon distribution is also in the RS case harder than the SM background. In the 
left panel of Fig. [7|we show da/dk±, having integrated over The steps occur when k± 
equals half the invariant mass of a RS-graviton, since this limits the photon momentum. 

The right panel of Fig. again shows bin-integrated cross sections, but now we have 
chosen 100 GeV bins around the first three resonances, a/I = m 1; m 2 , m 3 , since this is 
where most of the events will occur. We see a similar behaviour as in the ADD scenario. 
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Figure 7: Photon distributions for the RS model. Left panel: da/dk± (solid) with (from 
above) k/Mp\ = 0.1, 0.05, 0.01, mi = 1 and the SM background (dashed). Right panel: 
da/dk± (solid) with mi = 1 and k/Mj>\ = 0.05, integrated over 100 GeV bins around the 
first three resonances, with corresponding SM backgrounds. 
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6 Concluding remarks 



In summary, we have discussed photon Bremsstrahlung induced by the exchange of massive 
gravitons at the LHC Both the ADD and the RS scenarios have been considered. We 
found that three-body final states are likely to be detectable, and could be a valuable 
supplement to the two-body final states, for the purpose of detecting the effects from 
massive gravitons related to extra dimensions (if such exist). These configurations, of a 
hard photon associated with an electron-positron (or muon) pair in the opposite direction, 
should provide a striking signature at the LHC. Furthermore, the photon has a harder 
distribution in transverse momentum than the SM background. 

We have here focused on Bremsstrahlung at the LHC. At the Fermilab, the phenomeno- 
logy will be rather different. Because of the lower energy, and because of the different initial 
state, quark-ant iquark annihilation will be much more important, relative to the gluon- 
gluon initial state. Furthermore, there will in pp collisions be additional contributions to 
the forward-backward asymmetries, beyond those of the SM. Such asymmetries will be 
induced by interference between the C-odd exchange of a photon or a Z, and the C-even 
exchange of a graviton, as well as by the interference between graviton-exchange with 
initial-state radiation and graviton-exchange with final-state radiation. This effect, which 
is washed out at the LHC because of the symmetry of the initial state, is of course present 
also for e + e~ — ► // jSj, and will have an analogue in the three-body final states. We hope 
to discuss this effect elsewhere. 
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Appendix A: Convolution Integrals and Event Shapes 

In this Appendix, we define the convolution integrals and the event shapes used in our 
expressions. First, the basic convolution integrals are: 

I m i*)= f dyf g [\l-e*)f 9 [\l-^" 



s / \ V s 



Y 



s / \ V s 



W)= I dyf q [d-ey)fJJ-e-y\, (A.l) 



with Y = ~log(s/s), and where in the latter case, a specific quark flavor q is considered. 
The quark convolution integrals enter in the cross section weighted with different coupling 
constants and summed over quark flavors: 



Ibb(s) 


= 2^/ w -(s), 
q 




q 




= 2 

q 


JczCzis) 


q 


/d 7 d 7 (s) 


= 2j2QPqq(s), 

q 




= 2^2Q q v q I qq (s), 
q 


Id z d z (s) 


q 


J C 7 D Z (5) 


= IczD^s) = 2 Q 2 q a q I qq (s 
q 




= 2^2Q q v q a q I q q(s). 







(A.2) 

These integrals are labeled according to the sets of Feynman diagrams which are associ- 
ated with the couplings involved. Integrals involving SM couplings also enter in interference 
terms involving graviton exchange. Note the factor of two in the qq convolutions, which 
accounts for the fact that at a pp collider, either beam can provide the quark or the an- 
tiquark. All convolution integrals have been evaluated with CTEQ5 parton distribution 
functions [26 . 

The event shape distributions of the different contributions to the cross section can 
be expressed in terms of x±, x 2 and x% = 1 — x\ — x 2 - It is convenient to express these 
quantities in terms of the abbreviations: 

z a = 8xt - 12x 2 3 + 12x 3 ~ 3 > 

z b = 3x 2 3 (l - 2x 3 )(2x 2 3 - 2x 3 + 1), 

z c = 2x\( Xl + x 2 ) 2 (Ax 2 3 - 2x 3 + 1), 

^ = 6(1-2x3X4^-10x3 + 5), 

z e = 9x 2 3 (l - 2x 3 )(2x 2 3 - 6x 3 + 3), 

Zf = 8xj - 80x 3 3 + 180x l ~ UOx s + 35 > 
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Zg = + 2;r 3 - 1, 

z h = 2x\ - 6^3 + 3, 
Zi = 8x1 ~ ^ X 3 + 3, 
Zj = 8x1 — 10x 3 + 3, 

z k = 12x 2 3 ~ 8a; 3 + 3. (A.3) 

We first give the expression for initial-state radiation, expressed as an integral over 
cos#: 

f a (x 1 ,x 2 ) + a 2 (x 1 ,x 2 )cos 2 9 + a A (x 1 ,x 2 )cos i 9 
Xa{Xi > X2) = J d{C0Sd) ^(l-2x 3 )(l-cos^) ' (A ' 4) 



with 



a>o(xi, x 2 ) = -(x 1 - x 2 ) 4 z a - (x 1 - x 2 ) 2 z b + z c , 

a 2 (xi, x 2 ) = -(xi - x 2 ) 4 z d + (xi - x 2 ) 2 z e - x\z h , 

a 4 (xi, x 2 ) = (xi - x 2 ) 4 z f - (xi - x 2 ) 2 x\z d - x\z a . (A.5) 

For the initial-state SM background terms we have: 

Y f \ I A ( Co(x!,X 2 ) +C 2 (X 1 ,X 2 )C0S 2 ^ 

X c (x u x 2 ) = J d(cos9) 4{1 _ cos2e) , (A.6) 

with 

c (xi, x 2 ) = (xi - x 2 ) 2 z g + x\z h , 

c 2 (xi, x 2 ) = (xi - x 2 ) 2 z h + x 2 ^. (A.7) 

We impose a cut on the photon emission angle 9, given by \rj\ < 2.5, to render the integrals 
Xa and X c finite. 

For the remaining quantities, the integration over cos 9 is trivial, and one finds 
X B (x 1 ,x 2 )=X D (x 1 ,x 2 ) + S{Ai + 4) 



X D (xi,x 2 ) 



1-2X3 ' 

2(x 2 + x 2 ) 



(l-2x 1 )(l-2x 2 )' 



X AD (x u x 2 ) = (1 - 2x 3 ) 



(xi - x 2 ) 2 Zi + x\z 

4 

x 3 



3 



(xi - x 2 ) 2 z fc + x§(3 - 2x 3 ) 
^^'^) = x|(l - 2x 3 ) ' 

*CD0n,z 2 ) = ^±^. (A.8) 

Xo 
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Appendix B: Definition of y-cuts 



Here we shall define quantities where the xi,x% dependence (which determines event 
shapes) in the cross sections is integrated out. When we carry out these integrations, 
we will impose y-cuts. In the case of gluon-gluon fusion, we define 



These y-cuts will remove IR soft and collinear events where the photon has little energy, or 
its direction is close to that of an electron. The cut ys, which could be milder, will remove 
events where the two electrons are close. 

In the case of quark-antiquark annihilation, the approach is the same, so we will not 
write out the integrals here. However, for terms involving initial-state radiation and related 
interference terms which depend on S3 = (1 — 2xs)s, all factors involving S3 must be part 
of the integrand, since x\ and X2 are related to 23. 
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